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Abstract 

We derive a version of non-Abelian Stokes theorem for SU(2) gauge fields in which 
neither additional integration nor surface ordering are required. The path ordering 
is eliminated by introducing the instantaneous color orientation of the flux. We 
also derive the non-Abelian Stokes theorem on the lattice and discuss various terms 
contributing to the trace of the Wilson loop. 



Introduction 



The usual Abelian Stokes theorem relates the integral along closed curve C bounding some 
surface Sq and an integral defined on this surface. This version of the Stokes theorem is 
most relevant in physical applications since it allows to express the holonomy of gauge non- 
invariant electromagnetic potential via physically observable magnetic flux. The invention 
of non-Abelian gauge theories necessitated a non-Abelian generalization of the Stokes 
theorem. Nowadays, there are quite a few formulations of non-Abelian Stokes theorem 
(NAST) available (for review see, e.g., Ref. [I] and references therein). Generically there 
exist two principal approaches, the operator [2] — [7j and the path-integral jH] — [HI one. 

It is worth to mention that the central issue of any formulation of NAST is how 
to make sense of the path ordering prescription inherent to the non-Abelian holonomy 
(Wilson loop). In this respect neither operator nor path- integral approaches are helpful 
for concrete calculation of the Wilson loop. Indeed, while in the former case the path 
ordering is traded for much more complicated surface ordering prescription, in the latter 
case the additional path integral over auxiliary variables is introduced which cannot be 
calculated even approximately. 

In this paper we derive a new version of non-Abelian Stokes theorem focusing exclu- 
sively on SU(2) valued Wilson loops in the fundamental representation. Moreover, the 
central object of our discussion is the phase of the Wilson loop <p w . It determines the Wil- 
son loop trace, l/2TrVT = cos <p w , which is the only gauge invariant quantity associated 
with gauge holonomy. 

The basic idea is to introduce the instantaneous color orientation of the chromo- 
magnetic flux piercing the loop. Evidently this color orientation remains unknown until 
the Wilson loop is calculated by some other means. Nevertheless, it allows us to avoid 
the path ordering and represent the Wilson loop phase as an ordinary integral to which 
Abelian Stokes theorem applies. Furthermore, we relate the resulting surface integral 
with properties of gauge fields on this surface. 

As we noted above, in order to get explicitly the color orientation of the flux one has 
to calculate the Wilson loop first. In this respect our formulation is well suited for the 
lattice where the gauge holonomy is to be calculated numerically. We derive the lattice 
version of the non-Abelian Stokes theorem. Finally, we discuss the physical meaning and 
origin of various terms contributing to the trace of the Wilson loop and present the results 
of our qualitative numerical simulations. 

Non-Abelian Stokes theorem in the continuum limit 

Consider Wilson loop operator in the fundamental representation, W(T), evaluated on 
a closed contour C = {x^{t),t G [0;T],x /J (0) = x^(T)}, which is parameterized by dif- 
ferentiable functions x^{t). By definition the operator W(T) provides a solution to the 
first-order differential equation 

(V(t)| {idt + A) = 0, (1) 

(m\ = (m\w(t), (2) 



i 



W(t) = Texp{i f A{t)cLt}, (3) 
Jo 

where A = l/2a a A a ^[x{t)]x^(t) is the tangential component of the gauge potential, a a 
are the Pauli matrices, differential is taken with respect to the parameter t and ( ip \ 
is a vector in the spin-1/2 irreducible representation space (IRR) of SU{2) group. Since 
Eq. (JTJ) is nothing else but the time-dependent Schrodinger equation, the Wilson loop W(t) 
can be interpreted as a quantum mechanical evolution operator with the time-dependent 
Hamiltonian H = —A(t). Moreover, the corresponding state space coincides with spin- 
1/2 IRR, in which a convenient basis is provided by generalized (spin) coherent states ^3] 
(see, e.g., Ref. [TH] — [Ej for review). The spin coherent states {(n|} are parameterized 
by a set of unit three-dimensional vectors n, n 2 = 1 and in this basis arbitrary state ( ip I 
has a unique representation 

(^| = e iip (n| . (4) 
The action of SU(2) operator g reads 

(n\ g = e*e (H g \, (5) 

where the phase factor if depends on both n and g, ip = y)(g,n). Therefore, in the basis 
of coherent states Eq. (0) becomes 

(iP(t)\ = e lLp{t) ( n{t) | = (n{0)\W{t), (6) 

where without loss of generality we have taken y?(0) = or equivalently (^(O) | = (^(0) |. 

Eq. (^) imposes no restrictions on the initial vector ( n(0) |, which therefore can be 
taken arbitrarily. However, there exists a distinguished initial state which is of particular 
importance for the discussion below. Namely, let us take ( n(0) | to be the eigenstate of 
the full evolution operator 

(n(0)| W(T) = e^ (T) (n(0)|. (7) 

Note that generically l/2TrW(T) ^ ±1 and we assume this from now on. The gauge 
invariant trace of the Wilson loop is given by l/2Triy(T) = cosip(T). On the other 
hand, one gets from (I1I4|1 the following equation for the Wilson loop phase factor ^H] 

dt 



(f(T) = J (^(n \ A\n) — i{n\dt\n)^ 



Using standard properties of the spin coherent states [THj — [E| one can represent (jEJ) in 
vector-like notations 

" 1 f ^ . 1 



^ Tr W{T) = cos 
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[ nAdt + - [ n-ld^nx 3 v n\ d 2 a» v 
Jc 4 J Sc 



(9) 



where the vector field n(t) has been smoothly extended from the contour C into an ar- 
bitrary surface Sc bounded by C. Note that Eq. (jUJ) may be identically rewritten in the 
gauge invariant form 



]-TiW{T) = cos 



l -J^ {4, + n- [D^n x D v n] } d V 



(10) 
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where D"? — 5 



ab 



d u Al-d v Al-e abc A h „A 



o acb 74^ is the covariant derivative and - , .,, 

is the non-Abelian field strength. 

Let us emphasize that Eqs. (|9I1U|) cannot be used to actually calculate the Wilson loop 
since construction of the evolving states n{t) requires knowledge of the Wilson loop itself. 
Nevertheless, Eq. (jlOj) might be relevant for theoretical considerations since it represents 
the phase factor <p(T) as an integral of the field strength introduced first by 't Hooft [19] 
and Polyakov [20] in connection with monopoles. Eq. (fTT)|) is the non-Abelian Stokes 
theorem (NAST) to be discussed in more detail below. 

By construction the vector field n(t) is covariantly constant along the contour C, 
x^D^n = 0. Eq. © implies also that for any t £ [0;T] the state (n(t) | is an eigenstate 
of the Wilson loop calculated on C starting from the point In other words, (n(t) \ 

is an eigenstate of W + (t) W(T) W(t). The interpretation of n(i) is then straightforward, 
it is the instantaneous color orientation of the flux piercing the loop C. Since the contour 
considered is not infinitesimal the color direction of the flux is different at various points 
on C. 

The assignment of the vector field n(t) to a given closed path C is unfortunately not 
unique. The ambiguity comes from Eq. (JZj) which possesses two solutions with opposite 
sign of ft 



±n(0) | W(T) 



±n(0) 



However, this sign ambiguity is only global: if ( n(0) | is an eigenstate with ip(T) > 0, then 
for any t £ [0; T] ( n(t) | is an eigenstate of W + (t) W(T) W(t) with the same positive phase. 
Therefore there is only a global freedom to change n(t) — > —n(t) for all t simultaneously. 

Consider the infinitesimal closed contour SC X located in (//, v) plane at point x, which 
bounds the elementary surface element 6<t£ v . Since in this case W(T) = W(5C X ) = 
1 + ^<tF[ 1U 5<t x 11/ + o(5a) the eigenvalue problem (|7lllj) could easily be solved 



± 



F5a 
\F5a\ 



(p(5C x ) 



± \F5a\ 



± 



FSa) 



(12) 



where we have explicitly indicated the (fiu) dependence of the eigenvector and evident 
Lorenz indices have been suppressed. Therefore the infinitesimal version of non-Abelian 
Stokes theorem (jlUj) is (no summation over /i, v) 



-TrW{5C x ) 



cos 



In^F^Sa^ + o{5a) 



+ o(Sa 2 ), (13) 



where we have used the covariant constancy of n on 5C X . 

It is amusing to note that Eqs. (0 HUjl look similar to the non-Abelian Stokes theorem 
of Ref. |H] apart from the absence of the path integral over n(t) in Eq. (jlOj) . We conclude 
therefore that the path integral of Ref. [S] is exactly saturated by the two particular 
trajectories ±n(i). It is worth mentioning however that the construction of n(i), Eq. (|SjL 
makes no reference to the classical equations of motion and therefore Eqs. (|9I10|) do not 
correspond in general to any semi-classical approximation. 
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Notice that Eq. (fTUj) is still not uniquely denned. The point is that the vector field 
n(t) may be extended arbitrarily from C to Sc- The only requirement is that the extension 
C 3 n(t) — > n(a) G Sc must be continuous and the distribution n(a) must agree with 
n(t) at the boundary 5Sc = C. On the other hand, Eq. (fTUj) is applicable to any closed 
contour, in particular to any infinitesimal area element of Sc. Therefore, at every point 
x(a) G Sc we have a naturally defined direction n{a) which can be used to make Eq. ([10)1 
unambiguous. Then the only remaining problem is the choice of sign of n(a) since at every 
point x(a) G Sc considered separately there is no distinction between ±n(cr). Here the 
continuity requirement for the extension C 3 n(t) — > n(a) G Sc becomes crucial. Indeed, 
since the surface Sc is assumed to be regular (e.g., smooth and without self-intersections) 
the field F5a is continuous on Sc and therefore Eq. (JT2*Jl allows to define n(a) continuously 
as well (these arguments may fail for exceptional configurations which are not generic and 
which we don't consider for that reason). 

The sign ambiguity in the definition of n(a) is reminiscent to the model of Ref . [2JJ 122) 
(Alice electrodynamics). Indeed, the key feature of Alice electrodynamics is that the U(l) 
generator (n(a) in our case) is known only up to the sign. Moreover, the definition of the 
central object of the model, Alice string, is based on the continuity arguments similar to 
the above reasoning. However, the relevance of Alice electrodynamics to the 577(2) gauge 
theory is still unclear and we will not dwell on this issue. 

In the next Section we consider the non-Abelian Stokes theorem on the lattice. As a 
byproduct we also illustrate the appearance of various terms in Eq. fllOj) . 

Non-Abelian Stokes theorem on the lattice 

The derivation of lattice NAST begins from consideration of fundamental representation 
Wilson loop on the lattice 

N-l 

W = H Ui = U -U 1 -...-U N - 1 , (14) 

i=0 

where Ui G SU (2) are the link matrices parameterized as 

Ui = cosvi + i sinvi au i>+ = cos-Uj — i sin-Uj ff« ir , uf ± = 1 , (15) 

and the corresponding lattice path piercing the sites sq, si, sn-i is assumed to be closed 
and non self-intersecting. There is an analogous parameterization of the Wilson loop 

W = cos ip w + % sin (p w aw + = cos (p w — i sin cp w <jw_, w^. = 1 , (16) 

and therefore the state to be ascribed to the site so is ( w + \ = ( w + (s ) \ 

(w + (s )\W = e 1 ^ (w + (s )\. (17) 

Note that one could equally take the eigenstate (uL(s ) | instead (see discussion in the 
previous Section). As far as only a single Wilson loop is concerned, there are no much 
difference between the two choices and it is sufficient to take either of them. 
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Starting from (w + (so) | one constructs the corresponding eigenstates (w + (si) \, i — 
1, N — 1 in all other sites si, sn-i using Eqs. ()4|fjj) 

(v} + (si)\Ui = e**(u? + (s i+1 )|, i = 0,...,N-l. (18) 

Since the initial state was taken to be an eigenstate of W, the chain (jl8j) is closed, = s , 
(w + (sat) I = (w + (s ) | and one gets the following relation between the Wilson loop phase 
(f w and the phases 7i coming from links U{ [T8j 

<pw = 5^ 7i ■ ( 1Q ) 

We are in haste to add that, strictly speaking, Eq. (fT^j) is not entirely correct. The point 
is that the left hand side is always bounded, \<p w \ < vr while the sum on the right can take 
values outside the interval [— 7r;7r]. To be precise, Eq. (fTTIj) should express the equality 
of the phase factors e l<Pw , not the angles (p w by themselves. Therefore, there are terms 
2nk, k G Z missing in Eq. (|TT?|) . From now on in all the equations like ()19|) the mod 27T 
operation is always assumed and will not be indicated explicitly. 

In order to make one step further we need to consider in more detail Eq. (j!8|) . Consider 
an SU (2) operator U which upon acting on some initial state ( n\ \ brings it to another 
state ( n 2 \ 

(fii\U = e i7 (n 2 |, U = cos v + i sin v 3 u + , (20) 

where 7 is an additional phase which depends on both ( n\ \ and U. There are two 
eigenstates (u±\ of the operator U which form a complete basis in the spin-1/2 IRR 

(u±\U = e ±iv (u±\, 1 = \u + ){u + \ + (21) 

Using the resolution of unity (|2*T|) in (|2*U|) one gets 

e iu (ni|«+) (u+\ + e~ lv (n x \u-) (u-\ = e* 7 ( n 2 \ , (22) 

7 = v + il (n 1 ,u + ,n 2 ) = -v + £l (ni,u-,n 2 ) , (23) 

where Qo(fii, njv) is the oriented area of the spherical polygon on unit two-dimensional 
sphere S 2 with corners at the north pole of 5* 2 and at ni,...,n/v (in that order). In 
deriving Eq. (|2*3*|l standard properties of the spin coherent states have been used (see, 
e.g., Ref. [El— H3) 

(ni|n 2 )/|(n 1 |n 2 )| = e^C^.^). 

Note that in our normalization Area(S' 2 ) = 2n. In particular, there is no additional 1/2 
factor in front of Vt$. 

From very general arguments one expects that Eq. (|2()|) has an interpretation of ordi- 
nary rotation of the vector Hi around the axis u + by the angle v. However, as far as only 
initial and final states are taken into account the rotation operator remains in fact unde- 
termined. Indeed, there are infinitely many rotations which connect two given states. On 
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the other hand, among various SU(2) operators there is a distinguished unique G ni -> n2 
which describes the motion Hi — > H 2 along the shortest geodesic line connecting n%, H 2 

G ni ^ l2 = (Hi H 2 ) + 10 ■ [Hi x n 2 ] , ( n x \ G ni ^ n2 = e** ® 1 ™ (H 2 \. (24) 

The physical relevance of geodesic curves is widely known. The importance of the geodesic 
matrices G ni ^ n2 in the present context comes from the consideration of the following 
diagram 



u 



G, 



m 


— ► n 2 


ni 






G 






u 




mi 


— ► m 2 



Here the same operator U corresponds to the rotations Hi — > H 2 and mi — > m 2 . The 
diagram is closed by two geodesic matrices G mi ^ ni and G m2 ^ n2 . From the analysis of the 
diagram one obtains the relation 



G 



H U = UG 



m 2 ^n 2 ; 



(25) 



which follows from the fact that adjoint SU(2) action is equivalent to SO (3) rotation. 
Therefore the matrix UG m2 -> n2 U + is again geodesic with initial and final states being rhi 
and Hi and thus equals to G mi _ mi . Eq. allows to rewrite in the form 

Q (Hi,u ± , H 2 ) + £l (m 2 , u±, rhx) = Vt Q [m 2 , H 2 ) + Vt Q (Hi, m x ) . (26) 



Eqs. f!19l23l26|) are sufficient to derive the NAST on the lattice. We illustrate the 
derivation on the simplest example. Generalization to usual cubic geometry is straight- 
forward but technically notorious. Thus we will only discuss the final result. 





Fig. 1 

Consider the simplest non-trivial configuration of three links, Fig. 1, from which one 
could construct three different Wilson loops 



W = U Q U+ , Wi = U U+, W 2 = U X U. 



2 ) 



6 



W = W 1 W 2 . (27) 

Let tpi is the phase angle of the corresponding Wilson loop, 1/2 Tr Wi = cosifi. Applying 
the procedure of the previous Section to each Wilson loop separately, one gets six states, 
{n^ , Hg }, i — 0, 1, 2, sitting at the points A and B; pair of states with fixed i is assigned 
to the corresponding Wilson loops Wi. In particular, 

{n { 2 \ Wi = e^infl, i = 0,1,2. (28) 

Let us evaluate the phase angle (fo- According to ()19j23j) 

y? = v + J^^a > "b.+> ) - V2 + fi (^2 ) ,M 2 , + ,n^ ) ) , (29) 

where v% enters with minus sign because U 2 is conjugated in the definition of Wo- Using 
Eq. (J2fij) we can write 



(30) 



(fo = v + fio(^A } , f$) + flo(ny, u ,+, n% ] ) + Q {n [ g\ f$P) 

-v 2 + Qo{n ( g\n [ g } ) + fi (^2 ) ,^2,+,^A ) ) + ^o(^a\^a) 
The next step is to add zero in the form (see 



= Q (n A [ ) +^o(^2 ) ,^s ) ) + "1 + n (nA >"i,+>**b ) ~ v i + ^o("b\ Wi,+> ) ( 31 ) 
to Eq. (J30j) and collect various terms together using the relations like 

tpx = v + Qoin^, u ,+, fig) -vx + Qo{n { B, n^) , 

(32) 

0(A) = Qo(4 0) > + ^o(^ 1} , ^ 2) ) + CVn?, 4 0) ) , 

where f2(A) = fi(n^, n^, n^) is the area of spherical triangle constructed on the indi- 
cated three unit vectors at point A. In other words, Q(A) is just the oriented solid angle 
between the triple {n^ , , n^}. 

Therefore the final equation which relates the phase angles ifii, i = 0, 1, 2 is 

<p = (f! + <p 2 + V(A) + tt(B), (33) 

where il(A), Q(B) are the oriented solid angles between the triads {n^\ n^} and 
{fig , fig, fig}. Note the different ordering of states in Q(A) and £l(B), which corresponds 
to counting the outgoing flux at both points A and B. 

Let us emphasize that Eq. (}3*3*j) is valid irrespectively of the particular choice of the 
states {n^Ai^g} provided that the phases tpi are calculated according to (fT§j) or (J2Hj). 
It does not matter which particular solution of (|28p was taken to construct {w^ , n^} 
on each Wilson loop, Eq. (J3*3*j) remains formally the same with either choice. But this 
means that Eq. (J3l?|) is ambiguous because ipi changes sign when {n^,ng } are replaced 
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Fig. 2 



by {— nj^ , — fig }. In fact this is the same sign problem discussed previously. We will fix 
it after considering the continuum limit of Eq. (|33|) . 

The generalization of Eq. (f3l?j) to the case of usual cubic geometry is straightforward 
but technically involved. The final result 

p w = ^2 {px + ^ x + X] ax ( 34 ) 

xeSc x^Sc x&c 

is illustrated on Fig. 2. It is understood that phases ip w , ip x are calculated via Eq. 

is the phase of the large Wilson loop, l/2TrW(C) = cos<p w , where C is the planar 
2x2 closed contour, see Fig. 2, which bounds the surface Sc- The first term on the r.h.s. 
('dynamical part') is the sum of contributions coming from four 'internal' plaquettes 
belonging to Sc- In particular, 

-TrU Px = cos<p x , ^ = 0,1,2,3, 

where U Px is the corresponding plaquette matrix. The second term ('solid angle') comes 
from the points common to four different 'internal' plaquettes. We recall that application 
of NAST, Eq. (jHSJ), requires construction of four color vectors per plaquette situated at 
plaquette's corners. Therefore there are four unit vectors at the point x, see Fig. 2, and Q x 
is just the oriented solid angle between them. The third term ('perimeter contribution') 
is analogous to the second one. It accounts for the difference in color direction between 
the states on the nearest to the loop 'internal' plaquettes and the states on the loop itself. 
Technically a« is an oriented solid angle between the corresponding three vectors. 

Eq. has a simple physical interpretation. The magnitude of the total flux, ip w , 
piercing large closed contour C is the sum of a few terms. The first term sums up the 
magnitudes of elementary fluxes penetrating the surface Sc- Since the theory is non- 
Abelian each elementary flux has its own color orientation which is no less important than 
the flux magnitude (for flux piercing finite contour C the color direction of the flux varies 
along C). The other terms in Eq. (J34)) take into account the difference in color orientation 
of various fluxes on Sc as well as of the total flux piercing C. It is worth to mention that 
for pure Abelian fields (or for a SU (2) gauge copy of Abelian configurations) the second 
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and third contributions in (f3~4"j) vanish identically and one gets the usual Abelian Stokes 
theorem. 

Let us consider Eq. (|34j) in the limit of vanishing lattice spacing, a — > 0. The contri- 
bution of the first term was in fact already calculated in ()12I13|) 

i ^ i r - 

'dynamical part' = a 2 \_] -n x F /ll/ (x) + o(a 2 ) rs - / nF fMU d 2 a' J ' u , (35) 

Zr£. 2 4 JSa 



x<=S c 



where n x , x G Sc is given by (jl2j) . In order to get the continuum limit of the second term 
consider the point x G Sc and let (fAv) is a plane tangential to Sc at x. Then f2 x is the 
oriented solid angle between the four vectors 



n± = n x - a ■ (D^n x + D u n x ) , 

n 2 = n x - a- D v n x , 
n 3 = n Xl 

n A = n x - a- , 



(36) 



where n x is again given by |T2|) and D is the covariant derivative. It is straightforward 
then to evaluate Q x 

fl x = X -a 2 (n-[D^nx D u n}) + o(a 2 ) . (37) 

Therefore 

'solid angle' = ~ / n ■ [D^n x D u n] d 2 a^ . (38) 
4 he 

Unfortunately, there exists no simple expression for the third term, Eq. in the 

continuum limit. However, this is to be expected. Indeed, one can readily convince oneself 
that the meaning of the 'perimeter contribution' is to provide correct boundary conditions 
in Eq. (J33j) . In other words, the third term, Eq. pijl . guarantees that the vector field 
n(cr) G Sc agrees with n(t) G C on the boundary &Se = C. 

Combining Eqs. (j35l38j) one formally reproduces Eq. (jlOj) confirming that ()34j) is indeed 
the lattice formulation of the non- Abelian Stokes theorem ()10|). However, this conclusion 
relies heavily on Eq. (|37|) which is only valid if n X) x G Sc is continuous across the 
plaquette boundaries. This suggests a natural way to fix the relative sign of eigenstates 
on neighboring plaquettes analogously to the continuum considerations above. Namely, 
we propose to fix the particular distribution of eigenstates by the requirement that 

r = + Si^i ( 39 ) 

xdSc x&C 

takes the minimal possible value (it is assumed, of course, that eigenvectors at the bound- 
ary n(t) G C are held fixed from the very beginning). This prescription fixes completely 
and unambiguously all the states fz(cr) G Sc provided that the functional R has a unique 
minimum. The uniqueness of the minimum of R is a separate issue and we have no an- 
alytical methods to investigate it. However, at least numerically the minimum of (J3*9*|) 
might be approximated with high accuracy. 
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Numerical simulations 

In this Section we describe simple lattice experiments with Eq. (134)1 which we performed 
in pure SU(2) lattice gauge theory considered on 12 4 lattice at (3 = 2.4 using the standard 
Wilson action. 

Since the decomposition p4|) is gauge invariant (see discussion in previous Section) it 
is legitimate to ask what is the contribution of each term into the Wilson loop expectation 
value 

(exp{z^ }) ~ e~ TV ^, (40) 

where we have restricted ourselves to the consideration of rectangular T x R, T R 
loops only. Therefore the problem is to calculate 

(exp{ % <Px }) ~ e" TV ^ (R) , (41) 
xes c 

(exp{ i ^ V x }) ~ e - TV ° Md{R) , (42) 

(exp{ i^a x }) ~ e ~TV perim (R) _ ^ 

Notice that T, R dependence of the expectation values (j41H43|) is an ad hoc assumption 
which has to be checked separately. However, we have found that (j41M3j) indeed accu- 
rately describe numerical data. 

We calculated the expectation values (j40H43|) on 50 statistically independent config- 
urations using the spatial smearing algorithm (see, e.g., Ref. (22] f° r details). For each 
rectangular loop C = {T x R} Eqs. (I17ll8j) were applied to construct the eigenstates on C. 
The same procedure was used to build the eigenvectors {rip}, i = 0, 3 on each 'internal' 
plaquette p £ Sc (only surfaces with minimal area were considered). Finally, the func- 
tional (|39|) was minimized with respect to the inversions {n p 9} — > {—n p ^}, p £ Sq using 
a variant of the simulated annealing algorithm [21] and keeping the boundary conditions 
n £ C fixed. 

The results of our simulations are presented on Fig. 3, where circles represent the full 
heavy quark potential (jjO)), squares correspond the 'dynamical' part (j4T|) and finally dia- 
monds and triangles stand for 'solid' (}4*2*j) and 'perimeter' ()43|) contributions respectively. 
Note that the solid curves on Fig. 3 are drawn to guide the eye. 

There are few striking features of the expectation values (j41M3|) to be mentioned here. 
First of all, the 'perimeter' potential, Eq. ()43|). turns out to be practically /^-independent: 

V peHm (R) w const, (44) 

which might be an indication that the perimeter contribution drops out in the expectation 
value of the full Wilson loop (jjOj) . Secondly, both Vd yn {R) and V so i i( i(R) appear to be 
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linear at large distances, R>3, albeit with somewhat larger slope than the full potential 
V{R) « a SU ( 2 )R 

V dyn (R) » a dyn R, V soUd {R) = a soUd R, ° dyn & a "° l%d » 1.6. (45) 

&SU(2) C"St/(2) 

Although V dyn (R) deviates from the linear behavior at distances R < 3, K ^(i?) is 
rising strictly linear starting from the smallest possible distance R = 2. The existence 
of a linearly rising term in the heavy quark potential at short distances has been widely 
discussed in the literature, see, e.g., Refs. — [2H] and references therein. 

Finally, we emphasize that the expectation value of the full Wilson loop ()40j) is not 
factorizable into the terms (j41M3|) . It is clearly seen from Fig. 3 that 

V(R) ^ V dyn (R) + V solld (R) + V perim (R) + const (46) 

and therefore there are various interference terms contributing to V(R). The point how- 
ever is that the linear piece coming from V so ud might survive at small distances since the 
'solid angle' contribution ()38I42I45|) is formally not suppressed by the action even at very 
large (3. 



Conclusions 

We have derived a new version of the non-Abelian Stokes theorem for the Wilson loop in 
the fundamental representation of SU (2) gauge group. By considering the instantaneous 
color direction of the flux piercing the loop we were able to avoid the path ordering in the 
conventional definition of the Wilson loop operator. Moreover, this approach allows to 
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represent the phase angle of the Wilson loop as an ordinary one dimensional integral to 
which the usual Stokes theorem applies. Furthermore, we were able to relate the resulting 
surface integral with properties of non-Abelian gauge fields on that surface. 

Unfortunately, our formulation can hardly be called "theorem" because it does not 
help to calculate the Wilson loop itself. However, this drawback is not specific to this 
paper since other known variants of non-Abelian Stokes theorem are also not much useful 
for Wilson loop calculation. At the same time our construction is well suited for the 
numerical investigations. To achieve this goal we have also derived the non-Abelian Stokes 
theorem on the lattice and illustrated the origin and physical meaning of various terms 
contributing to the trace of the Wilson loop. 
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